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Abstract. The influence of non-uniform distribution of nuclei on crystallization 
kinetics of amorphous materials is investigated. This case cannot be described by 
the well-known Johnson-Mehl-Avrami (JMA) equation, which is only valid under 
the assumption of a spatially homogeneous nucleation probability. The results of 
computer simulations of crystallization kinetics with nuclei distributed according to 
a cluster and a hardcore distribution are compared with JMA kinetics. The effects 
of the different distributions on the so-called Avrami exponent n are shown. 
Furthermore, we calculate the small-angle scattering curves of the simulated 
structures which can be used to distinguish experimentally between the three 
nucleation models under consideration. 



PACS numbers: 07.05.T, 61.12.E, 61.43, 61.43.B, 05.40 
1. Introduction 

The properties of metallic glasses and other amorphous materials may be impaired 
by even small amounts of crystalline phases. Crystallization can also improve the 
properties of some amorphous materials, e.g. glass ceramics. So the understanding 
of the crystallization process is very important. The analysis of experimental data is 
often made within the framework of the Avrami theory |Q, |^ by means of the Johnson- 
Mehl-Avrami equation, which gives a relation between the fraction of transformed 
(i.e. crystallized) material x(t) and the time at constant temperature. An equivalent 
approach was made by Kolmogorov |^ . 
Both models are based on the following assumptions: 

(i) Crystallization is considered in an unlimited medium. 

(ii) Nucleation of crystals begins at time t = and occurs in a non-crystallized region. 
The nucleation rate per unit volume, a{t), is assumed to be independent of the 
coordinates. 

(iii) The growth of crystals ceases at the points of mutual impingements, whereas it 
continues unchanged elsewhere. Before they touch each other, the crystals have 
a geometrically similar and convex shape (the Avrami approach is restricted to 
spherical crystals). 

(iv) The growth rate v{t) in a given direction is the same for all crystals and depends 
only on time. 



2 



Based on these assumptions, one can derive the following exact relation: 



X{t) = 1 - exp 



where 



V{T,t)=Vo 



-J a{T)V{T,t)dT 




v{t')dt' 



(1) 



(2) 



Vo is a form factor, in the case of spherical crystals Vq = 47r/3. An integration 
of equation ([l]) is only possible by making specific assumptions about the time 
dependence of the nucleation rate a(t). 

If both V and a are independent of time (interface controlled growth and continuous 
nucleation), then 



X(t) = 1 - exp 



(3) 



If only the growth rate v is constant and all crystals are formed simultaneously at i = 
with a mean number density /? (instantaneous nucleation), the resulting nucleation rate 
a{t) — (35 (t) can be substituted into equation (0). In this case 



X{t) 1 - exp [-fiVovH^ 



(4) 



Avrami proposed that for a three-dimensional nucleation and growth process with 
constant or decreasing nucleation rate, the general relation 



X{t) - 1 - exp(-fct") 



(5) 



should describe the crystallized volume fraction. Equation (||) is the so-called Avrami 
equation with the Avrami exponent n, where i < n < A. 

This equation is often used to analyse experimental data by means of a logarithmic 
plot, where In {— ln[l — x(i)]} is plotted versus h\[t). The slope of the resulting 
straight line is the Avrami exponent n, which describes crystallization kinetics. But if 
the previously mentioned assumptions are not exactly satisfied, the resulting Avrami 
exponents n may be misleading. We study cases where some of the assumptions (|-0) 
are violated and no analytical results for the crystallization kinetics are available. 
If the chemical composition of the two phases involved in the transformation is 
different, the growth rate is diffusion controlled. In this case, the radii r of spherical 
crystals grow according to 



r{t',t) = g ■ Vt-t' 



(6) 



where g is a constant and t and t' are the observation and nucleation times, 
respectively. Substituting (^ into equation (|l]) leads to an Avrami-exponent n with 
1.5 < n < 2.5. But if continuous nucleation occurs, assumption (|i^) is violated, as 
the growth rate decreases with increasing life time of the crystal. As a consequence, 
crystals nucleated at different times have different growth rates. In the Johnson-Mehl 
method, the growth law (^ allows "phantom crystals" (crystals that are nucleated in 
an already crystallized area) to outgrow the real crystal. Therefore, the crystallized 
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volume fraction is overestimated. 

The assumption (|i|) has no physical reasons. In practice, nucleation may occur 
preferentially in certain macroscopic regions (e.g. grain boundaries), so that the 
nucleation probability becomes dependent on the coordinates 

Furthermore, the grains have to reach a critical size rcrit to start growing. This critical 
radius is determined by the free energy and the surface energy. Therefore, in a sphere 
of radius rcrit only one crystal can nucleate This also violates assumption (^). 
In the present paper the influence of non-uniformly distributed nuclei and critical 
radii on crystallization kinetics is studied by computer simulation. We use methods of 
stochastic geometry ^ 0], namely the germ-grain models. The models are explained 
in section ^ and the simulation technique is pointed out in section ^ The results are 
presented in section 0. 



2. The models 



A germ- grain model is defined by a point field Pi, P2, ■•• with density A and a series of 
grains Ai, A2, ... with finite size. The complete model A is formed by the union of the 
grains A„ shifted to the points P„. Due to the diversity of point fields and types of 
possible grains there is a great variety of germ-grain models. In our study, we always 
use spherical grains of radius r^. 

To model the pure JMA case with assumptions (|-[^ fulfilled, we use a Poisson model 
for the underlying point field. The two fundamental properties of this point field are: 

• The number, N{G), of points lying in an arbitrary region G with volume V{G) is 
a random variable. The probability P of finding n points in the region G is given 
by the Poisson distribution 

P{N{G) ^n)= [^^(f )]" exp[-Al^(G')] n=l,2,... . (7) 

• Considering disconnected regions Gi,G2,..., the numbers N{Gi), N{G2), ■■■ are 
independent random variables. 

Germ-grain models with an underlying Poisson point field and overlapping grains are 
also called Boolean models. For these models it is possible to calculate the volume 
fraction xa analytically ||, ||, 0, |^ : 

XA = 1 - exp[~XViA)] . (8) 

Here, A is the number density of the Poisson point field and V{A) is the mean volume 
of the grains. Equation (||) is equivalent to equations (|| - ^) if the corresponding 
nucleation and growth laws are inserted. 

To model a system with increased nucleation probability in certain regions, we use 
a cluster point field. The nuclei are distributed uniformly and independently within 
spheres of radius R^. These spheres are distributed according to a Poisson point field 
with parameter Apar, the numbers of points within the spheres are Poisson distributed 
with mean value A'ci- Hence, the density of the cluster point field is Ad = Apar-^ci- 
An appropriate set of parameters for characterizing the cluster model is {Xc\, Nc\,c) 
where c = 2Rc\/fi. Here, fi is the mean distance of the midpoints of neighbouring 
clusters given by [Q: 

ri = (1^) r(4/3)A-i/' « 0.554A-y3 . (9) 
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For c 3> 1 the point field approaches a Poisson point field of density Ad- 
A hard-core ■point field is used to force a certain minimum distance between the nuclei. 
The distance between any points of the model with density Ahc is forbidden to be 
smaller than a given value i?hc- The essential property of the structure is described 
by the packing fraction phc — 47r/3-R^^Ahc- For phc — > the point field approaches a 
Poisson point field. 

Neither for the cluster model nor for the hard-core model analytical expressions are 
known that describe the volume fractions. 

The small-angle scattering intensity I{q) per unit volume is given by 

oo 

I{q) =4^1 r2[CA(r) - ^iWl^^^dr . (10) 



To calculate the small-angle scattering intensities of the germ-grain models, the 
covariance CA(r) of the systems is needed. This is the probability P of two random 
points ri , r2 with distance r both lying in the region covered by the model: 

CA{r) = P{ri e A, f2 e A) , r =1 ri - r2 I . (11) 

For Boolean models an analytical expression Q exists: 

Ca = 2xA - 1 + (1 - XAf exp [A7^r)] (12) 

where 7°(?') is the mean distance probability function averaged over all spheres with 
density f{x) of the radii distribution: 

oo 
r/2 

The constructional details of the point fields mentioned above are explained in the 
next section. 



3. The simulation technique 

The nucleation and growth processes are simulated in a cube of unit lenght Lq and 
volume Vo = Lq. All lengths are scaled to Lq. To model an infinite structure, periodic 
boundary conditions are applied. We use spheres of equal (instantaneous nucleation) 
or different (continuous nucleation) size as grains. They grow according to the specified 
growth law and nuclei are generated randomly according to the underlying point field. 
In the case of instantaneous nucleation (INST), all grains start growing at i = 0. If 
continuous nucleation (CONT) is considered, in every evolutionary step dt a mean 
number of nuclei starts growing according to the nucleation rate a{t). In this case, 
the nuclei that are created in an already transformed area have to be ommitted. 
After every time step, the volume fraction XA of the system is calculated. Optionally, 
the covariance Ca (r) of the structure can be calculated at a given volume fraction. 
In every simulation, 500 - 700 nuclei are generated. To limit the influence of statistical 
fluctuations, the whole procedure is repeated 10 - 40 times and the average of the 
relevant quantities is evaluated. 
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3.1. Construction of the point fields 

In a first step of the simulation tlic underlying point field has to be created. 
To generate a Poisson point field, the number N^^ of nuclei is drawn from a Poisson 
distribution with parameter AVq . Then A^nu points with cartesian coordinates equally 
distributed in (0,Lo) are created. 

A cluster point field is build by a two-step procedure. First, a Poisson point field of 
parent points with parameter Apar is created in Vq. In a second step, spheres of radius 
Rci are attached to these parent points. In each of these spheres, another Poisson 
process with parameter Nd is created. Omitting the parent points, the remaining 
nuclei obey a cluster distribution with density Ad = Apar^ci- 

To construct a hardcore point field, single points with equally distributed coordinates 
are generated subsequently. New points are accepted only if their distance to all exist- 
ing points is greater than 2i?hc- This procedure is repeteated until the desired number 
of nuclei (Poisson deviated with parameter XhcVo) is reached. 



3.2. Nucleation and growth 

After the definition of the point fields, the time evolution of the system starts. 

In the case INST, all predefined grains start growing at t = 0. In every time step, 

their radii are calculated as r(t) = w(r)dr. All grains have the same size. 

If continuous nucleation is simulated, in every time step iVact nuclei start growing. iVact 

is Poisson deviated with the parameter a(t)dt. The radii of the grains are different 

now, depending on the life time ti of each individual grain i: 

ti 

n{t.i) = / v{t)At . (14) 





3.3. Calculation of the volume fraction and the covariance 

To calculate the volume fraction, a fine grid of Attest test points is constructed in the 
unit cube. The coordinates of these test points are determined using a quasi-random 
sequence according to Sobol ||l^, |ll|, |l^. The volume fraction is given by the number 
of test points A^in that lie within the area covered by the spheres: 

XA = ^ • (15) 

To evaluate the covariance, the distances between all the test points are calculated. 
The numbers N(r) of test points that lie within discrete distance intervalls between r 
and r + dr are determined. Then, the covariance is given as 

where Ni^(r) denotes the number of test points in the distance intervall (r, r -I- dr) 
that are covered by a sphere. 

The small-angle scattering intensities can now be calculated according to equation 
( p^ ) . To solve the integral, we use a Fourier transfom method and fit a polynomial to 
the discrete values of CA(r) — Xa order to obtain the function to be integrated. 
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3.4- Test of accuracy 

Due to the finite and discrete nature of the simulation, two sources of systematic errors 
have to be considered: 

• The size of the system (i.e. the number of grains) has to be large enough in order 
to describe an infinite system. 

• In the simulation of CONT, the nucleation proceeds in small but finite time steps 
dt. 

In order to check the accuracy of the simulation, we performed calculations using the 
Poisson point field and compared the results for the volume fraction and the covariance 
with the exact equations (||) and (Q), respectively. 

In figure |l|, the results of simulations in the case of instantaneous and continuous 
nucleation are shown. The simulations were performed 10 times with a step size 
dt = 10 (INST) resp. dt = 1 (CONT)| and 5 • 10^ test points were used. An Avrami 
analysis by means of linear regression yielded an Avrami exponent n = 2.96 for the 
simulation of instantaneous nucleation, which is in good agreement with the exact 
value n — 3. 

In the case of continuous nucleation, the value of the Avrami exponent is n = 3.95, 
compared with the exact value n — A. 

The simulated small-angle scattering intensities of instantaneous nucleation 
(parameters as above) at different volume fractions are shown in figure |2[ Here, 
quantitative differences between the simulated values and the exact ones calculated 
according to (^, (|2|) and (^ occur, although the covariance values are in quite 
good agreement, see figure ^. In both cases, the same numerical integration method 
was used. Because of the multiplication of Ca(?') — Xa by r in equation ([lO|), very 
small differences between the simulated covariance values and the exact ones at large 
r-values yield substantial differences after integration. 

Hence, with the present accuracy of our method, only qualitative statements 
concerning the small-angle scattering curves are possible. But the main features of 
the scattering intensities at different volume fractions are represented properly. The 
curve with low volume fraction shows well-resolved maxima and minima as the single 
spheres are still nearly isolated. With increasing volume fraction, the amplitudes of 
the oscillations decrease, at xa = 0.9 there are only weak ripples left. 

4. Results and discussion 

We investigated the dependence of crystallization kinetics on the spatial distribution of 
the nuclei in the case of continuous nucleation (CONT) and instantaneous nucleation 
(INST). Additionally, for CONT the case of diffusion controlled growth was surveyed. 
The small-angle scattering intensities were calculated for instantaneous nucleation in 
order to check if it is possible to distinguish between the several grain distributions 
by means of small-angle scattering methods. 

Instantaneous nucleation 

In the case of instantaneous nucleation, we checked the infiuence of cluster- and 
hardcore- model on crystallization kinetics. The results were compared with the 

t The time scaling does not have any influence on the resulting Avrami exponent. 
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Figure 1. Comparison of the simulated volume fractions (•) with the exact JMA- 
equation {-) in the case of instantaneous nucleation (INST) and homogeneous 
nucleation (CONT). The error bars denote the standard deviation of the results. 
INST: A = 500, V = 0.002, CONT: a = 5, v = 0.002, dt = 1. 



corresponding simulated values for the Poisson model. We used constant growth rates 
V for calculations^. 

Furthermore, the covariance values and the resulting small-angle scattering intensities 
were evaluated. 



Results on crystallization kinetics In figure^, Avrami plots of a cluster model 
and a hardcore model are compared with that of a Poisson model with equal point 
density A. It is clearly shown that the distribution of the nuclei according to a 
cluster model leads to a reduced Avrami exponent n compared with the value n = 3 
expected for the Poisson model. For small values of c, the simulated values deviate 
from a straight line in the Avrami plot, which means that the crystallization kinetics 
cannot be represented by an exponential law according to in these cases. If c is 
small enough, n raises again and the deviations from the exponential law decrease 
again. The values of the simulated Avrami exponents (drawn from a linear regression 
analysis) are shown in table |l| for two simulations with (Ad = 50, Nd = 10, c) and 



I Note that in the case of instantaneous nucleation time dependent growth rates can be reduced to 
this case as all grains start growing at the same time. See M . 
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Figure 2. Comparison of the simulated small-angle scattering values (dashed lines) 
with exact results (solid lines) for Poisson distributed nuclei at different volume 
fra<;tions. The parameter a denotes the radius of the grains. \ = 500. 

(Aci = 25, TVci = 20, c), respectively. Simulation results in the case of a hardcore 

Table 1. Avrami exponents for several parameter values of liie cluster distribution 
and the hardcore model in case INST with A = 500. Poisson model; n = 2.97 



Cluster, N^i = 10 Cluster, Af^i = 20 Hardcore 
c n c n phc ^ 



2.66 


2.36 


2.11 


2.08 


0.26 


3.56 


2.00 


2.26 


1.58 


2.01 


0.13 


3.27 


1.33 


2.26 


1.06 


2.07 


0.06 


3.07 


0.67 


2.46 


0.53 


2.33 


0.02 


2.99 


0.33 


2.43 










0.17 


2.64 










0.08 


2.80 










0.04 


2.88 











model with Ahc = 500 and several hardcore radii i?hc are also listed. The simulation 

of the corresponding JMA case yields an Avrami exponent n = 2.97. According to 
these simulation results, a distribution of the nuclei according to a hardcore model 
leads to an increased Avrami exponent n compared with the one expected by the JMA 
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Figure 3. C omp arison of tlie simulated covariance values (o) with exact results (— ) 
according to (|l2[). Parameters as in figure ^ (INST). 

theory. A deviation from the hnear behaviour in the Avrami plot can also be observed 
for large values of phc- 



4- 1-2. Results on small-angle scattering data To check the possibility of 
distinguishing experimentally between the distributions under consideration, we 
investigated the small-angle scattering curves at several volume fractions. 
The scattering curves of the Poisson model were already shown in figure |^. Figure ^ 
shows scattering curves of a germ-grain model with underlying hardcore distribution. 
For low volume fractions, the curves exhibit a significant first peak. This peak is 
characteristic for a hard-sphere model with non-overlapping spheres (see, e.g. 0). 
With increasing volume fraction this peak disappears as the structure of the system 
is now far away from the structure of the generating nuclei and the overlapping of 
the grains becomes larger. The scattering curves of a cluster model with different 
volume fractions are shown in figure |^. Here, at low volume fractions no sharp peaks 
are present. With increasing volume fraction, the amplitude of the oscillations first 
increases and then decreases again. 

Considering these results, it should be possible to distinguish between the three 
distributions of the nuclei by using small-angle scattering and observing the whole 
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Figure 4. Influence of the parameter c of the cluster model and the parameter p^c 
of the hardcore model on the Avrami exponent in case INST. The straight lines are 
linear regression values. A = 500, A^cl = 10. 

crystallization process. 

4-2. Continuous nucleation 

The same investigations as in the case of instantaneous nucleation were also made 
for continuous nucleation and linear growth rate v = const.. In the case of diffusion 
controlled growth, we analysed the error that is made by applying the JMA equation. 
In both cases, predefined point fields of density A = 750 and a nucelation rate a = 5/dt 
were used for the simulation. 

In figure 0, the influence of a cluster and a hardcore point distribution on crystallization 
kinetics is shown (the given parameters describing the point fields refer to the 
predefined nuclei). Further simulation results are listed in table |^. As in the case 
of INST, a nuclei distribution according to the cluster model results in a decreased 
Avrami exponent compared with the JMA case and in deviations from the exponential 
behaviour. A hardcore distribution leads to an increased Avrami exponent, but the 
differences to the Poisson model are not as distinct as in case INST. 
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Figure 5. Simulated small-angle scattering curve of a harcore model with different 
volume fractions. A^c = 500, -Rhc = 0.045. The parameter a denotes the radius of 
the grains. 



Table 2. Avrami exponents for several parameter values of the cluster distribution 
and the hardcore model in case CONT with a = 5/dt. Poisson model: n = 3.95 



Cluster, = 10 


Cluster, N^i = 20 


Hardcore 


c n 


c 


n 


Phc n 


4.57 3.72 


3.63 


3.56 


0.27 4.04 


3.05 3.59 


2.42 


3.33 


0.20 3.99 


1.52 3.40 


1.21 


3.14 





4-.3. Discussion of crystallization kinetics 

The deviations from the Poisson model in the cases INST and CONT can be explained 
by looking at the derivation of the JMA equation (see, e.g. To calculate the trans- 
formed volume fraction, a so-called extended volume is introduced, which is simply the 
sum of the volumes of all grains without considering their mutual overlappings. To 
get a connection between this extended volume and the real transformed volume, as- 
sumption ([^) is applied. 

If the grains are distributed according to the cluster model, their overlapping is under- 
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Figure 6. Simulated small-angle scattering curve of a cluster model with different 
volume fractions. Ad = 50, A^cl = 10, c = 2.0. The parameter a denotes the radius 
of the grains. 



estimated by the JMA model and hence the transformed volume fraction is overesti- 
mated. For c » 1 (clusters lose their cluster-like nature) , n is close to the JMA- value. 
With decreasing c of the underlying cluster point field, the Avrami exponent n gets 
smaller and substantial deviations from the exponential behaviour of the JMA-kinetics 
occur. If c — > 0, the Avrami exponent increases again and the Avrami plot shows a 
linear behaviour. In this case, the clusters are widely spaced and act like single 
Poisson-distributed nuclei. In between these two limiting cases (c ^ 1 and c ^ 0), 
deviations from the JMA-behaviour occur. Values of c with approximately 1 < c < 2 
give the maximum deviation (minimum n). 

On the other side, a hardcore distribution leads to a smaller overlapping compared 
with a uniform distribution. The transformed volume fraction is larger than in the 
pure JMA case, since more of the space nuclei grow into is empty. Our simulations 
show an increase of the Avrami exponent n with increasing packing fraction phc of the 
underlying point field. On the other hand, in the case of phc n reaches the value 
of the pure JMA-case. Unfortunately, with our present algorithm we could not reach 
the limiting case of a close-packing of the underlying point field. 
As outlined in section [11 in the case of diffusion controlled growth the crystallized 
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Figure 7. Influence of the parameter c of the cluster model and the parameter p^c 
of a hardcore model on the Avrami exponent in case CONT. The straight lines are 
linear regression values. A = 750, N^i = 20, a = 5/dt. 

volume fraction is overestimated. To check for this, we performed simulations with 
an underlying Poisson distribution of the nuclei. In a first series, the phantom nuclcii 
that nucleated in an already crystallized region were discarded. Afterwards, they were 
treated like regular grains and contributed to the volume fraction. Doing so, we could 
estimate the error that is made in applying the JMA equation on diffusion controlled 
growth. The calculation of the differences Aabs of the volume fractions yielded differ- 
ences Aabs < 9 • 10""^. These results are in good agreement with simulations made by 
Shepilov and Bochkarev [^ . 



5. Conclusions 

Our simulations concernig the dependence of crystallization kinetics on the spatial 
nuclei distribution clearly showed that the analysis of experimental data by the JMA 
equation must be done with care. If the nuclei are not distributed equally, the use 
of the JMA equation can yield substantially wrong results. Therefore it should be 
checked if the JMA equation is applicable. One possibility to do so is the use of small- 
angle scattering. The scattering curves of the investigated nuclei distributions differ 
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clearly from one another, especially in an early stage of the crystallization process. 

On the other hand, the simulations showed that the error that is made by applying 
the JMA equation on diffusion controlled growth processes with continuous nucleation 
can be neglected. 
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